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A contribution  of  COIL  Laboratory  in  Prague  to 
the  Chemical  Oxygen-Iodine  Laser  research  and  development 

Jarmila  Kodymova*  and  Otomar  Spalek 
Department  of  Gas  Lasers,  Institute  of  Physics  of  Academy  of  Sciences 

Prague,  Czech  Republic 

ABSTRACT 

The  key  results  gathered  in  the  COIL  Laboratory  of  the  Institute  of  Physics  AS  in  the  Czech  Republic  since  1985  to  date 
on  the  experimental  and  theoretical  investigation  of  Chemical  Oxygen-Iodine  Laser  (COIL),  and  related  problems  are 
reviewed  in  a certain  context  of  historical  perspective  of  the  COIL  research  and  development. 

Keywords:  chemical  lasers,  oxygen-iodine  laser,  singlet  oxygen,  atomic  iodine 


1.  INTRODUCTION 

Laser  emission  at  1.315  pm  from  magnetic-dipole  transition  between  spin-orbit  levels  of  the  ground  state  configuration 
of  iodine  atom 

iePm)(F=i)^l{2Pil2HF=4)  (1) 

was  first  observed  during  flash  photolysis  of  alkyliodides,  and  since  its  discovery  by  Kasper  and  Pimentel1  in  1964,  the 
iodine  photodissociation  laser  has  been  developed  into  systems  capable  of  very  high  energies.  The  Prague  Asterix  Laser 
System  (PALS)  is  currently  one  of  the  largest  such  facilities  in  the  world,  delivering  on  the  target  a focud  power  density 
up  to  1016  W/cm2  in  400  ps  (at  the  fundamental  wavelength).2  Lasing  mechanism  of  the  Chemical  Oxygen-Iodine  Laser 
(COIL)  follows  the  same  physical  principle,  only  the  2PU1  upper  level  is  pumped  via  the  energy  of  singlet  delta  state  of 
molecular  oxygen,  02(1Ag).  Unique  physical  and  chemical  properties  of  this  metastable  species3  that  is  able  to  transfer 
its  energy  to  iodine  atom  in  a near-resonant  process,  and  can  be  generated  chemically  with  very  high  efficiency  have 
brought  about  the  exceptional  position  of  COIL  among  other  chemical  lasers.  The  pumping  mechanism  discovered  by 
Derwent  et  a/.4  in  the  seventies,  is  represented  by  the  reverse  reaction  process  with  temperature  dependent  equilibrium 
constant 


I(2P3/2)  + OAg)  I(2Pi/2)  + 02(3Ig)  Kc<r  kjk~  0.75  exp(401.4/T)  (2) 

Twenty-five  years  have  passed  in  2002  since  McDermott  and  his  colleagues5  from  the  Air  Force  Weapons  Laboratory 
(the  Air  Force  Research  Laboratory  nowadays)  demonstrated  the  first  COIL  generation  on  the  miliwatt  level.  A 
multihundred-kilowatt  version  of  the  COIL  facility  has  been  developed  and  tested  for  the  Airborne  Laser  program6 
since  that  time.  In  the  meantime,  a tremendous  theoretical  and  experimental  work  has  been  devoted  to  understand 
comprehensively  this  laser  system  to  bring  it  to  the  present  technological  level.  The  Prague  laboratory  joined  this 
investigation  in  1985.  Theoretical  and  experimental  efforts  have  been  focused  on  partial  problems  in  the  COIL 
operation  that  have  concerned,  for  example,  the  C^Ag)  generation  for  both  subsonic  and  supersonic  COIL,  the 
chemical  generation  of  atomic  iodine  for  a COIL,  laser  performance  characteristics  in  CW  and  pulsed  regime,  etc.  A 
brief  cross-section  through  results  of  this  investigation  is  presented  below. 
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2.  SINGLET  OXYGEN  GENERATION 


2.1.  Mechanism  of  O^Ag)  loss  in  a low-pressure  chemical  generator 

The  25th  commemorative  COIL  session  permits  mentioning  the  interesting  result  following  from  our  “pioneer”  work 
done  on  the  02(1Ag)  generation  by  the  reaction7 

Cl2  + HO/  ->  O^Ag)  + H+  + 2C1'  fe  = 10n  cm3morls*1  , (3) 

and  the  mechanism  of  02(1A)  loss  in  a low-pressure  bubbler-type  (sparger)  generators8-10  employed  for  driving  the 
COIL  in  subsonic  flow  regime.  The  02(1Ag)  formed  in  the  liquid  film  surrounding  the  bubbled  chlorine  through  the 
BHP  has  to  diffuse  into  the  gas  bubble  during  the  time  that  the  bubble  is  moving  up  through  the  head  of  liquid  (basic 
solution  of  hydrogen  peroxide,  BHP).  The  gas  pressure  (consisting  of  reacting  chlorine  and  formed  singlet  oxygen) 
within  the  two-phase  (g-1)  bubbling  mixture  in  the  reactor  column  is  substantially  increased  due  to  a hydrostatic 
pressure  of  the  liquid  column  than  in  the  gas  space  above  the  liquid  head.  For  example,  if  the  02  pressure  in  the 
generator  plenum  is  1 Torr,  then  the  8 cm  liquid  head  provides  the  hydrostatic  pressure  of  8 Torr  on  gas  bubbles  at  the 
reactor  bottom.  By  this  effect,  the  loss  of  02(1Ag)  in  gas  bubbles  at  the  bottom  is  64times  higher.  A mathematical  model 
for  estimation  of  these  loss-processes  in  sparger-type  generators  was  developed  considering  a variation  in  the 
volumetric  fraction  of  gas  in  g-1  mixture  with  a distance  from  reactor  bottom,  as  well  as  the  effect  of  water  vapour 
pressure  in  gas  bubbles  on  the  02(!Ag)  loss.9  This  simplified  model  is  based  on  the  mass  balance  for  02(*Ag)  over  a 
differential  layer  dy  of  the  g-1  mixture  following  from  the  relation 

-NAndrj  = r4pSdy  (4) 

Na  - the  Avogadro  constant,  n - molar  flow  rate  of  chlorine  or  oxygen,  rj  - 02(lAg)  yield  from  layer  dy  (y  - vertical 
coordinate  specifying  the  distance  from  reactor  bottom),  S - reactor  cross  section,  r4  - rate  constant  of  reaction  (5). 11 

02(lAg)  + 02(lAg)  02(,Sg)  + 02(3Eg)  ks  = (2.7±0.4)xl0-17  cmWV1  (5) 

The  p factor  is  the  volumetric  fraction  of  gas  in  the  g-1  mixture  evaluated  from  the  height  of  two-phase  (foamy)  mixture 
(h)  and  a liquid  height  without  gas  in  the  reactor  (ho),  P - 1 - hjh.  There  are  two  main  02(*Ag)  deactivation  processes 
in  these  generators.  First,  the  02(*Ag)  in  the  bubble  diffuses  back  to  the  liquid  surrounding  it,  is  attached  and 
deactivated.  Secondly,  a deactivation/pooling  02(1Ag)  self-quenching  process  occurs  in  the  gas  phase  (bubble)  itself  by 
the  reaction  (5).  This  process  is  pressure  dependent  (quadratic  in  02(]Ag)  concentration),  and  consequently  strongly 
accelerated  within  g-1  mixture  due  to  the  hydrostatic  pressure  of  liquid  (in  the  above  given  example,  by  two  orders  of 
magnitude).  Contrariwise,  if  the  bubble  residence  time  in  BHP  is  too  short,  not  all  Cl2  is  converted  into  02(1Ag)  and 
generator  yield  is  also  low.  The  results  of  this  research  are  remembered  in  Fig.  1. 


Fig.  1.  02(1Ag)  yield  as  a function  of  BHP  height  in  reactor.  1-  experimental  curve, 

2 - calculated  curve  including  02(!  A)  loss  by  reaction  (5),  3 - calculated  curve  including 
another  02(lA)  losses  (in  liquid  phase,  quenching  on  reactor  wall,  transport  in  duct) 
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It  is  evident  that  this  type  of  singlet  oxygen  generator  can  operate  at  moderate/low  pressure,  and  only  way  to  achieve 
high  mass  flow  is  to  increase  the  generator  surface  area  resulting  in  a larger  02  collection  plenum.  This  finding  was 
afterwards  considered  in  designs  of  sparer-type  generators. 

2.2.  Generation  of  Oa^Ag)  in  a high-pressure  jet  SOG 

Development  of  the  supersonic  COIL  during  early  ninetieth  called  for  a high-pressure  02(!Ag)  generators.  We 
concentrated  on  the  investigation  of  the  jet-type  SOG  that  was  set  up  first  in  the  Samara  COIL  laboratory.12  A large 
surface  area  of  liquid  in  the  form  of  liquid  jets  for  performing  the  reaction  (3)  allowed  a substantial  reduction  of  the 
generator  volume  (dimensions),  and  consequently  the  02(1Ag)  loss  in  the  gas  phase.  Velocity  of  liquid  jets  is  very  high 
(5-20  m/s),  which  ensures  a fast  renovation  of  the  gas-liquid  boundary,  and  so  minimizes  H02'  ions  depletion  from  the 
surface  layer.  It  contributes  to  lowering  the  02(!Ag)  loss  also  in  the  liquid  phase.  Our  effort  in  the  investigation  of  this 
generator  was  focused  on  some  problems  in  operation  of  this  generator.13'1^  Experimental  results  obtained  on  the  small- 
scale  pilot  device  concerned,  for  example,  the  influence  of  different  configurations  and  diameters  of  holes  (0.3, 0.5,  and 
0.8  mm)  in  the  BHP  injector  on  the  jets  stability,  and  a detailed  investigation  of  water  vapour  estimation  in  gas  exiting 
the  generator.  We  searched  for  conditions,  at  which  the  water  vapour  content  was  minimized. 

Water  formed  during  02(1Ag)  generation  causes  together  with  peroxide  vapor  a significant  problem  in  the  COIL 
operation  because  it  quenches  the  excited  state  of  atomic  iodine  in  very  fast  processes17 


H20  + I(2P|/2>  ->  H20  + I(2/*3q) 

ki,  - (2.3±0.3)xl0'12  cmV 

(6) 

H202  + I(2P,,2)->H202  + I(2/>3/2) 

kn  = (1.0±0.2)xl0'n  cmV1 

(7) 

In  a supersonic  laser,  water  can  cause  also  undesirable  fluid-dynamic  effects  due  to  its  condensation  during  the  gas 
expansion  cooling.18  Water  vapor  concentration  in  gas  was  estimated  by  the  relation  (8)  derived  from  the  mass  balance 
of  02(lZg)  state.19 

CH20  = *4  CA2/  k9  Cx  (8) 


This  species  is  formed  in  the  pooling  reaction  (5),  and  is  very  efficiently  quenched  by  water  molecule.20 

O 2(lSg)  + H20  0Z('AJ  + H20  k 9 = 6.7xl0"12  cmV1  (9) 

Water  partial  pressure  was  evaluated  from  the  emissions  of  O^’lg)  and  02(1Ag)  states  for  different  BHP  composition, 
BHP  jets  temperature,  varying  Cl2  flow  rate,  and  different  BHP  injector  configurations.  The  experimentally  obtained 
data  were  compared  with  calculated  values  of  the  saturated  water  vapour  pressure  above  the  BHP  (for  given  BHP 
compositions  and  temperatures).  The  both  values  were  nearly  the  same.  The  relative  water  vapour  content  was  well 
decreasing  on  increasing  the  generator  pressure.  It  was  the  experimental  proof  that  employing  the  high-pressure  jet  SOG 
in  the  COIL  operation  reduces  a detrimental  effect  of  water. 

2*3.  The  Einstein  coefficient  for  02(1Ag)  fundamental  emission 

Work  on  the  development  of  high-pressure  jet-type  chemical  generator  of  02(lAg)  and  the  above-described  parametric 
investigations,  and  simultaneously  in  the  literature  opened  discussions  on  the  right  value  of  the  Einstein  coefficient  for 
02(*Ag)  emission  at  1.27  pm  (,4 -coefficient)  triggered  our  interest  in  verification  of  this  coefficient.  Mlynczak  and 
Nesbir1  challenged  the  Badger  et  aV s22  value  of  2.58  x 10*4  s'1  when  they  reported  a new  value  of  1.47  x Hf4  s'1.  The 
radiative  lifetime  of  02(*Ag),  X&9*,  -65  min  (a  half-lifetime  -45  min)  considered  so  far  in  the  COIL,  increased  to  -1 13 
min.  A detailed  insight  into  this  problem  revealed  that  published  values  of  the  Einstein  >4 -coefficient  evaluated  by 
means  of  various  theoretical  approaches  and  experimental  techniques  varied  in  a wide  range  resulting  in  xArad  from  -53 
to  -151  min.23  24  This  fact  could  make  the  evaluation  of  02(*Ag)  concentration  in  COIL  questionable. 

To  verify  a rightness  or  faultiness  of  the  ,4-value  used  so  far,  the  way  of  employing  two  independent  experimental 
techniques  simultaneously  was  chosen.25  The  electron  paramagnetic  resonance  (EPR)  spectroscopy  was  one  method, 
and  the  optical  emission  spectroscopy  was  the  second  one  supplemented  with  the  mathematical  model  developed  for  the 
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evaluation  of  02(!Ag)  concentration  from  the  emission  measurements.  In  this  model,  we  assumed  that  02(lAg)  molecules 
are  distributed  uniformly  in  the  optical  cell  emitting  photons  in  all  directions  with  a rate,  r , 

r = A&c&s  = A&P&  / RT  (10) 

- Einstein  coefficient  for  spontaneous  emission  of  02(lAg),  PA-  O^Ag)  partial  pressure. 

The  Oa^Ag)  concentration,  c^,  was  evaluated  from  the  measured  signal  of  the  fundamental  emission  at  1.27  pm  and 
declared  sensitivity  of  the  photodiode,  and  the  calculated  ratio  of  the  02(1  Ag)  emission  detected  by  the  photodiode  and 

the  emission  of  all  02(1Ag)  molecules  present  in  the  detection  cell.  The  proposed  mathematical  model  was  used  to 
evaluated  this  ratio,  in  which  a photocurrent  of  the  photodiode  was  directly  proportional  to  a value  of  the  Einstein  A - 

coefficient  substituted  into  relation  (10).25  A quantitative  evaluation  of  the  02(1Ag)  concentration  by  means  of  the  EPR 
spectroscopy  was  based  on  the  “double  integral  method”  of  EPR  spectra  processing  by  the  relation25 

oo  B 

caEPR=  Jfcjfj  y*  dBdB  (11) 

-00-00 

By  comparison  of  the  02(*Ag)  concentrations  obtained  contemporary  by  both  experimental  methods  and  supplemented 
by  the  mathematical  model,  a new  value  of  the  Einstein  ^-coefficient  of  (2.24  ± 0.4)  x W4  s'1  was  obtained.  The 
corresponding  02(lAg)  radiative  lifetime  is  - 74  min.  Shortly  after  our  first  reporting  this  result,  Newman  et  al 26 
published  the  value  of  (2.19  ± 0.07)  x 10*4  s'1  giving  the  rAad  ~ 76  min. 

2.4.  Generation  of  02(1AS)  in  a discharge  plasma 

The  theoreticians  of  our  department  performed  during  the  eightieth  the  fundamental  investigation  on  the  physical 
kinetics  of  high-frequency  (hf)  glow  discharge  in  oxygen  with  respect  to  generation  of  the  a*Ag  state  for  pumping  the 
iodine  laser.22  Coming  out  from  the  properties  of  the  electron  gas  described  by  solution  of  the  Boltzmann  equation  for 
the  electron  distribution  function,  the  02(lAg)  concentration  in  the  plasma  was  calculated  with  precondition  that 
population  of  this  electronic  level  is  driven  by  electron  collisions  only.  The  upper  concentration  limit  of  02(*Ag) 
achievable  in  this  reaction  system  was  evaluated  from  the  statistical  balance  between  the  production  and  loss  processes. 
The  calculated  relative  population  of  a'Ag  state,  XJN,  is  shown  in  Fig.  2 as  a function  of  the  hf  field  parameters  E^fN 
and  co  (E^-  effective  value  of  hf  field,  N - concentration  of  neutral  molecules,  co  - frequency).  It  can  bee  seen  that  the 
yield  of  02(*Ag)  in  hf  discharge  plasma  in  pure  oxygen  can  hardly  exceed  15%  in  the  chosen  range  of  parameters. 
Moreover,  it  is  the  maximum  02(*Ag)  concentration  that  might  be  still  reduced  by  collisions  with  negative  ions  present 
in  the  plasma,  e.g.27 


02(*Ag)  + O'  ->  03  + e k\2~  3xlO‘10  cm3  s*1  (12) 

Outside  the  region  of  active  discharge,  the  relative  02(*Ag)  concentration  is  even  lower  (only  3 to  10%  in  the  pressure 
range  of  1-10  Torr).  According  to  this  studies,  the  hf  discharge  unlikely  might  produce  the  a*Ag  state  in  the  amount 
sufficient  to  sustain  the  operation  of  iodine  laser  pumped  by  the  energy  of  such  generated  singlet  oxygen.  This 
conclusion  was  supported  by  a subsequently  performed  analysis  by  Swift28  of  the  kinetic  processes  that  accompany  the 
02(!Ag)  generation  in  microwave  discharge  plasma  in  pressure  range  of  1-10  Torr.  By  these  calculations,  a maximum 
attainable  relative  concentration  of  02(!Ag)  was  10%. 

A current  experimental  research  on  plasma  discharge  generation  of  02(*Ag)  for  pumping  the  iodine  laser  (Discharge 
Oxygen-Iodine  Laser,  DOIL)  is  performed  by  our  colleague  J.  Schmiedberger  in  his  long-time  co-operation  with 
Japanese  Institutes.  The  last  results  were  obtained  with  the  generation  of  radio-frequency  plasma  of  02+N2+N0  mixture 
in  the  hollow  electrode  jet  generator.29  The  generator  is  schematically  shown  in  Fig.  3.  The  plasma  jet  was  extra  chilled 
reactively  with  transversely  injected  mixture  of  Ar+N02.  Adding  NO  to  oxygen  enhanced  the  02(1Ag)  production  up  to 
20%  through  the  energy  transfer  from  electronically  excited  state  of  NO  molecule  to  02.  Another  increase  was  achieved 
when  N02  was  injected  into  the  reaction  system,  which  caused  via  certain  kinetic  processes  a lowering  in  a content  of  O 
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and  03  species  in  the  gas  mixture.  The  02(1Ag)  with  32%  yield  was  generated  at  pressure  of  0.43  Torn  In  combination  of 
this  rf  plasma  generator  with  a transversally  injected  molecular  iodine,  atomic  iodine  emission  at  1.315  pm  was 
detected  at  the  laser  output  mirror  by  the  spectrum  analyser.  This  emission  was  recalculated  to  the  output  power  of 
~3  nW.29 


TO  LASER 


Fig.  2.  Calculated  relative  concentration  of  the  02(1Ag)  as  Fig.  3.  Radio-frequency  plasma  jet  generator  of  C^Ag) 
a function  of  E^N  for  several  values  of  the  parameter  co/N : 

1 -2.65xl0'8;  2 “4.7x1  O'8;  3-9.4xl0*8;  4“1.33xl0'7; 

5 - 1.88x1  O'7;  6 - 2.65x1  O'7  [cmV]. 


3.  ATOMIC  IODINE  GENERATION 

A generation  of  atomic  iodine  in  the  ground  state,  I(2p3/2)»  directly  in  the  COIL  medium  to  replace  its  conventional 
production  from  dissociated  molecular  iodine  by  02(!A)  opens  up  a way  for  potential  improvement  of  the  COIL 
performance.  The  first  is  the  removal  of  the  02(*A),  energy  to  dissociate  I2,  allowing  this  energy  to  be  extracted  as 
power.  The  second  is  the  decrease  in  the  molecular  weight  of  the  injectant,  allowing  for  faster  diffusive  mixing,  and 
avoiding  negative  effect  of  I2  on  the  I(2Pi/2)  quenching.  Also  the  technical  problems  induced  by  special  treatment  of  I2 
vapour  could  be  diminished. 

With  these  needs  in  mind,  our  efforts  has  been  focused  currently  on  the  gas  phase  chemical  method  of  atomic  iodine 
generation  employing  a principal  reaction 

X + HI  — » I + HX  X = F or  Cl  (13) 

The  reaction  (13)  represents  in  fact  a multiple  gas  phase  chemical  processes  that  occur  in  the  overall  reaction  system.  It 
involves  initially  the  generation  of  F and  Cl  atoms,  respectively,  in  the  reaction  of  molecular  fluorine,  F2,  or  chlorine 
dioxide,  C102,  with  nitrogen  oxide,  NO.  Additionally,  each  of  these  processes  is  composed  of  multiple  individual 
chemical  reactions  describing  production  and  loss  processes,  including  species  contained  in  the  COIL  medium.  And 
finally,  this  set  of  processes  occurs  within  a flowing  gas  with  the  exothermicity  of  the  reactions,  which  brings  into  the 
system  a complication  with  coupling  the  chemical  kinetics  to  the  fluid  dynamics. 

The  goals  of  our  work  on  this  problem  are:  1/  to  provide  an  overview  of  the  current  theoretical  understanding  of  the  two 
complex  reaction  systems,  and  to  fit  them  to  experimental  conditions  characteristic  for  a COIL,  2/  to  test  experimentally 
the  both  reaction  systems  in  a flow  reactor  (a  pilot  small-scale  device)  to  prove  their  suitability  for  a realistic  COIL,  and 
3/  to  test  prospectively  the  developed  method  of  atomic  iodine  generation  directly  on  an  adapted  supersonic  COIL 
device  built  in  our  laboratory.  Work  towards  the  first  goal  has  involved,  to  date,  a model  based  on  the  generalized  one- 
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dimensional  (1-D)  flow  development,  exploring  the  chemical  processes  within  the  reaction  systems.  The  calculations 
have  been  performed  for  physical  conditions  (pressure,  temperature,  and  flow  rates)  applicable  to  a subsonic  region 
(i.e.,  upstream  of  the  nozzle  throat)  of  the  lkW-class  supersonic  COIL  operated  in  our  laboratory.  The  effects  of  initial 
concentration  ratio  of  reactants  and  different  configuration  of  reactants  injection  were  considered.  A rate  of  I atoms 
generation  has  been  investigated  in  a flow  of  nitrogen  (as  a buffer  gas),  and  a reactive  stream  containing  02(lA).  A 
reaction  mechanism  for  the  kinetics  of  I atoms  generation  in  the  (C1+HI+02(1A))  and  (F+HI+02(1A)),  respectively, 
system  has  been  developed  and  sensitivities  to  key  processes  in  this  mechanism  have  been  considered.  The  modeling 
provided  the  concentration  profiles  of  all  relevant  reaction  components  and  the  temperature  profile  along  the  reaction  path. 
The  results  of  this  relatively  rough  estimation  showed  that  I atoms  can  be  generated  via  F atoms  with  the  yield  of  up  to 
70%,  however,  on  a rather  long  reaction  path  (20-40  cm).  The  reaction  system  via  Cl  atoms  provided  by  these 
calculations  I atom  yield  80%  on  substantially  shorter  reaction  path  (-4  cm).  For  this  reason,  the  later  system  would 
meet  better  the  COIL  mixing  conditions.  Atomic  iodine  could  be  generated  near  the  critical  nozzle  plane  of  the 
supersonic  COIL  or  during  the  transonic  mixing.  The  chemical  kinetics  of  fluorine  system  is  however  simpler,  and  all 
reactants  are  commercially  available  in  gas  cylinders.  A detailed  description  of  results  on  the  1-D  modeling  of  both 
reaction  paths  for  chemical  generation  of  atomic  iodine,  and  their  interpretation  were  partly  published30*31,  and  some 
advanced  results  can  be  found  in  the  paper  of  this  Proceedings.32  A more  sophisticated  estimation  of  optimal  conditions 
for  I atoms  generation  needs  the  2-D  or  3-D  modeling,  involving  also  the  fluid  dynamics  effects.  This  work  is  planned 
for  a near  future. 

Work  towards  the  second  goal  above  has  involved  to  date  the  experimental  investigation  of  I atoms  generation  via  Cl 
reaction  system.  A design  of  the  small-scale  flow  reactor,  and  the  experimental  conditions  were  developed  on  the  basis 
of  the  results  of  the  above-mentioned  modeling.  Preliminary  experiments  provided  atomic  iodine  with  the  molar  flow 
rate  up  to  160  pmol/s,  and  I yield  over  70%  (in  nitrogen).33  Effects  of  initial  ratio  of  reactants,  and  the  way  of  their 
mixing  were  studied  during  the  following  on  investigation.  The  yields  of  atomic  iodine  in  the  range  of  70-100%  (related 
to  HI  or  C102  flow  rates)  were  attained.  Some  advanced  results  can  be  found  in  the  paper  of  this  Proceedings.34  The 
experiments  on  atomic  iodine  generation  in  the  02(lA)  stream  are  under  way. 


4.  REPETITIVELY  PULSED  SUBSONIC  COIL 

A period  of  the  early  ninetieth  signified  one  of  the  meaningful  millstones  in  the  COIL  R&D  - a demonstration  of  ability 
to  operate  a COIL  in  the  pulsed  regime.  A prestige  of  COIL  among  other  laser  systems,  particularly  technological  ones, 
so  has  increased.  The  magnetic  gain  modulation  in  a COIL  could  make  possible,  e.g.,  a more  efficient  non-linear 
conversion  of  the  wavelength  1.315  pm  to  the  second  and  higher  harmonics.  A slow  (non-pulsed)  control  of  3-4  gain 
could  be  utilized  for  a continuous  regulation  of  CW  power  (within  0-100%  of  the  maximum  value)  without  any  change 
of  flow  or  other  operation  parameters.  Further,  it  could  be  employed  as  a method  of  the  output  power  stabilization  by 
means  of  the  negative  feedback  controlling  the  magnetic  field  intensity  in  the  gain  region. 

Besides  previously  investigated  pulsed  methods  applicable  to  the  COIL,  like  a pulsed  UV  photolysis  of  alkyl  iodides  in 
gas  containing  singlet  oxygen* 8 5,  and  electric  discharge  initiation  of  I atoms  production  from  alkyl  iodides6'8*36,  a Q- 
switching  technique  either  mechanical37  or  magneto-optical38  was  also  applied  to  a COIL.  The  optical  switching  utilizes 
the  phenomenon  of  suppressing  the  gain  on  the  basic  of  the  Zeeman  effect.39  Owing  to  a low  gain,  low  output  coupling, 
and  near  threshold  generation  in  a COIL,  the  0-switching  showed  a disadvantage  due  to  an  increase  in  intracavity 
losses. 

The  pulsed  technique  investigated  in  our  laboratory  is  also  based  on  principle  of  the  magnetic  suppression  (modulation) 
of  gain,  but  the  Zeeman  effect  is  induced  by  externally  applied  magnetic  field  on  the  gain  region.40  The  small  signal 
gain  on  the  strongest  spectral  line  3-4  observed  in  the  laser  generation  in  zero  magnetic  filed  is  weakened  in  growing 
magnetic  field,  and  due  to  above-mentioned  physical  characteristics  of  COIL,  rather  low  magnetic  fields  are  needed  to 
quench  the  laser  generation  at  1.315  pm.  Such  a threshold  or  under-threshold  condition  can  be  expressed  by  the 
relation41 


(2gUB)L)/(S0+S')<\ 
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gu(B)  - small-signal  gain  at  magnetic  field  B,  L - gain  length,  S0  - intracavity  losses,  Sc  - resonator  output  coupling. 
Peak  power  enhancement  over  cw  power  can  be  achieved  in  the  pulsed  mode  when  pulsed  power  extraction  times  are 
short  compared  to  the  time  that  it  takes  to  fill  the  outcoupled  gain  volume  with  fresh  “frier  (in  a subsonic  system 
usually  '-lO'3  s).  Energy  of  the  single  pulse  is  given  by  the  average  CW  power  (with  zero  magnetic  field),  and  the 
optimal  repetition  rate.  A pulsed  length  and  its  enhancement  over  CW  power  depend  also  on  a switching  off  rate  of  the 
magnetic  field.  The  gain  volume,  from  which  power  can  be  extracted,  is  limited  by  deactivation  rate  of  both  the  singlet 
oxygen  and  excited  atomic  iodine.  Provided  that  this  rate  is  fixed  for  a particular  temperature  in  the  cavity  region,  the 
way  to  increase  this  volume  and  thus  increase  the  peak  power  enhancement  is  to  increase  the  velocity  of  the  flowing 
laser  gases,  and  so  to  increase  the  concentrations  of  02(1Ag)  or  atomic  iodine  or  both.  Experimental  verification  of  the 
proposed  pulsed  technique  was  performed  on  the  small-scale  subsonic  COIL  device  generating  the  CW  output  power  up 
to  several  tenth  watts.  It  was  modified  to  apply  magnetic  field  to  the  gain  cavity  from  two  rectangular  coils  attached  to 
the  laser  body  made  of  Plexiglas  (see  Fig.  4). 
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Fig.  4.  Scheme  of  set  up  for  magnetic  pulsed  regime  of  COIL 


A typical  oscillogram  of  the  periodically  pulsed  COIL  is  shown  in  Fig.  5 that  was  recorded  for  15  W of  the  average  CW 
power.  The  peak  enhancement  corresponded  to  the  power  of  48  W,  i.e.,  the  enhancement  ratio  was  3.2.  A length  of 
single  pulse  was  ~200  ps  at  the  rate  of  magnetic  filed  pulsing  of  0.5  G/ps.41  A maximum  peak  power  enhancement 
achieved  on  this  experimental  configuration  was  ~3.5  only  with  the  repetition  rate  of  ~1  kHz.  It  was  given  by  the 
relatively  long  switching  off  time  of  the  used  magnetic  field  generator.  Consequently,  extraction  of  the  pulse  power 
could  not  be  fully  developed. 


Fig.  5.  Oscillogram  of  the  repetitively  pulsed  COIL, 
upper  trace  - pulsed  magnetic  filed  (1  ms/div,  peak  current  250  A); 
lower  trace  - laser  pulsed  power  (CW  power  15  W,  pulsed  power  48  W) 
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Concurrently  with  our  work,  the  magnetic  field  nullified-switched  gain  technique  was  developed  at  the  AFRL  by  Hager 
et  al42  This  technique  works  on  principle  of  the  gain  suppressing  below  the  lasing  threshold  by  a static  magnetic  field 
that  is  followed  by  applying  a fast  rising  current  pulse  to  cancel  the  cavity  magnetic  field.  A switching  off  the  magnetic 
field  was  very  fast  in  this  way  (-0.5  ps)  resulting  in  the  laser  pulse  extraction  with  duration  less  than  20  ps,  and  a peak- 
CW  power  enhancement  ratio  of  12.7. 

Experimental  data  on  estimation  of  the  “3-4  gain  quenching”  magnetic  field  intensity,  called  a threshold  magnetic  field, 
were  also  collected  during  this  investigation.  It  is  defined  as  g{B(h)  = 0).43  The  magnetic  field  was  continuously 

increased  up  to  level  when  the  output  laser  power  quenched.  The  values  of  g(Bth)  evaluated  in  dependence  on  iodine 
molar  flow  rate  in  the  cavity,  and  for  several  values  of  output  coupling  parameter  are  plotted  in  Fig.  6.  This  dependence 
showed  a linear  character  in  the  overall  tested  parameters,  even  though  the  dependence  of  laser  power,  as  well  as  output 
coupling,  on  iodine  molar  flow  rate,  showed  generally  maximum  for  a certain  value  of  iodine  concentration  in  the 
cavity. 


Fig.  6.  Dependence  of  threshold  magnetic  field  on  iodine  molar  flow  rate  for  different  output  couplings 

It  can  be  concluded  from  these  results  that  a certain  value  of  threshold  magnetic  field  is  affected  by  a set  of  physical  and 
operation  parameters  together,  particularly  by  iodine  molar  flow  rates,  small  signal  gain,  threshold  gain,  and  output 
couplings.  In  spite  of  the  fact  that  a mutual  relation  among  these  parameters  is  rather  complex,  we  proposed  using  a 
dimensionless  parameter,  the  normalized  gain  defined  as  the  ratio  of  threshold  and  small-signal  gain  (g^/g),  to  evaluate 
from  the  experimental  data  the  threshold  magnetic  field  for  given  conditions.  The  threshold  gain  includes  the  effect  of 
output  coupling,  intra-resonator  losses,  and  gain  length.  Iodine  flow  rate,  Doppler  and  pressure  broadening  of  the  laser 
transition,  and  also  mixing  efficiency  affect  the  small  signal  gain.  The  dependence  of  normalized  gain  on  the  threshold 
magnetic  field  is  shown  in  Fig  7. 


Fig.  7.  Normalized  gain  ( gth/g0 ) in  relation  to  threshold  magnetic  field 
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5.  OPERATION  OF  1 kW-CLASS  SUPERSONIC  COIL 


5.1.  Performance  characteristics 

A small-scale  5-cm  gain  length  supersonic  chemical  oxygen  iodine  laser  driven  by  a jet  singlet  oxygen  generator  (jet 
SOG)  with  a closed-loop  liquid  flow  system  has  been  developed  since  1996,  and  laser  characteristics  have  been 
investigated.  Parameters  of  the  principal  device  parts  are  similar  to  the  VertiCOIL  developed  in  the  AFRL  except  the 
rotating  disks  generator.  The  PC  on-line  with  designed  hardware  and  software  nearly  fully  automatically  controls  the  all 
system  and  records  32  essential  laser  parameters.  An  iodine  vapour  in  helium  is  injected  into  the  primary 
oxygen/helium  flow  in  a subsonic  region  5 to  7 mm  upstream  from  the  nozzle  throat  plane.  A single  horizontal  slit 
configuration  of  nozzle  is  used  for  the  adiabatic  expansion  of  gas  flow.  A distance  between  the  sonic  plane  and  optical 
axis  of  the  laser  resonator  is  adjusted  from  35  to  55  mm.  A stable  configuration  of  the  resonator  with  5 cm-diameter 
inner  mirrors  is  employed.  The  mirrors  are  protected  by  helium  atmosphere.  The  nozzle  ramps  and  mirrors  holders  limit 
the  multimode  output  beam  size  to  3.7  cm  in  the  flow  direction  by  1.5  cm  between  the  ramps. 

Few  experimental  results  are  shown  to  represent  state  of  the  art  of  the  investigation  performed  on  this  device  to  date. 
The  jet  SOG  is  operated  typically  at  pressure  7.5  to  9.0  kPa  (60  to  70  torr),  with  helium  admixed  to  chlorine  in  the  ratio 
mostly  1:2.  A pressure  in  the  generator,  and  consequently  02(1Ag)  yield  could  be  easily  controlled  by  throttling  the  gas 
at  the  generator  exit.  It  can  be  seen  from  one  example  of  the  generator  performance  record  shown  in  Fig.  8.  An  increase 
in  the  yield  with  pressure  decreasing  results  from  a shorter  residence  time  of  gas  in  the  generator,  and  also  a lower 
Oii'Ag)  partial  pressure.  Both  effects  diminish  02(  'Ag)  loss  in  the  self-quenching/pooling  reactions11 

O^'Ag)  + ©.(’A,)  ->  02(3Zg)  + O^Sg)  *15=2.7xlO  '7cmV  (15) 

02('Ag)  + 02('Ag)  ->•  02(3Zg)  + O/Eg)  ki6  = 1.7xl0'17  cm3  s1  (16) 


Fig.  8.  Pressure  in  jet  generator  and  02(lAg)  yield  controlled  by  gas 
throttling  at  generator  exit.  Position  of  throttle  valve  flap  in  tens  of  degrees 
is  denoted  by  TV  numbers  (TVO  - closed,  TV9  - opened  entirely). 

Cl2  flow  rate  of  40  mmol/s,  Hepnm  flow  rate  of  40  mmol  /s. 

The  jet  SOG  driving  the  supersonic  COIL  provided  typically  rather  high  02(1Ag)  yields  (see  Table  1). 

Tab.  1 


PEen.kPa 

7 

8 

9 

10 

02(‘ac)  yield,  % 

80 

77 

70 

63 

94 
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Gas  dynamic  conditions  in  the  whole  laser  system  were  tested  during  a „cold  flow  run“,  i.e.,  with  gaseous  nitrogen 
instead  of  chlorine  was  introduced  into  the  primary  flow.  The  average  Mach  number,  Mi,  for  a subsonic  flow  ahead  of 
the  sonic  throat,  and  the  supersonic  flow  in  the  resonator  region,  M2,  were  evaluated  from  pressure  (Pi ,2)  in  the 
respective  regions  of  the  device,  gas  molar  flow  rates,  (w),  physical  properties  of  the  flowing  gas  (molecular  weight, 
MW , and  adiabatic  constant,  k),  stagnation  temperature,  and  the  flow  cross  section.44  The  local  Mach  number,  M2\  in 
the  supersonic  region  was  evaluated  from  the  pressure  ratio  PSta/Ppit  (Pstat  - static  pressure  measured  at  the  cavity  wall, 
Ppit  - stagnation  pressure  measured  by  the  Pitot  tube  in  the  cavity  centreline  placed  55  mm  from  sonic  throat) 45  The 
results  are  summarized  in  Table  2. 

Tab.  2 


^N2» 

mmol/s 

«He, 

mmol/s 

Pi, 

Pa 

Pi, 

Pa 

^ Pit, 

Pa 

MW, 

kg/mol 

K 

M, 

m2 

M* 

21 

90 

2055 

174 

1457 

8. 6x1  O’3 

1.57 

0.41 

1.94 

237 

Laser  output  power  was  investigated  at  various  iodine  flow  rates  and  output  couplings.  One  example  of  such 
dependence  is  shown  in  Fig.  9.  A simplified  Rigrod  theory46  for  saturation  intensity  was  employed  to  estimate  the  small 
signal  gain  and  saturation  intensity  from  these  experimental  data.  The  gain  of  1.29  x 10'2  cm'1,  and  the  saturation 
intensity  of  5.4  kW/cm2  were  evaluated. 

The  effect  of  titration  ratio  [l2]/[02(1Ag)]  on  laser  power  was  searched  also  for  different  output  couplings.  Fig.  10 
illustrates  it  for  two  values  of  generator  pressure.  The  highest  power  was  mostly  achieved  for  [I2]/[02(IAg)]  = (5±1)%. 
The  optimum  value  of  1.7%  obtained  on  the  VertiCOIL  system47  was  significantly  lower,  however  this  discrepancy 
diminishes  to  3.1%  when  the  titration  ratio  is  recalculated  by  dividing  the  total  [02]  in  VertiCOIL  with  02(1Ag)  yield 
(0.55). 


r^rrmoi/s 


Fig.  9.  Laser  power  as  a function  of  iodine  flow  rate. 
Output  coupling  2.83%  (Ti+T2),  40  mmol  Cl2/s, 

80  mmol  Heprim/s,  40  mmol  Hescc/s. 


Fig.  10.  Laser  power  as  a function  of  titration  ratio 
[I2]/[  C^Ag)]  at  two  pressures  in  generator 


The  optimum  dilution  ratio,  [Hesec]/[I2],  at  maximum  attainable  laser  power  was  found  in  the  range  from  1:30  to  1:70. 
These  measurements  were  directed  to  get  efficient  mixing  of  the  primary  gas  flow,  02(1Ag)+HePriin,  and  the  secondary 
gas  flow,  I2+Hesec.  It  was  estimated  by  means  of  the  relative  penetration  factor,  tc^-  rin The  penetration  parameters, 
;rand  flfoii,  describing  the  mixing  both  gas  flows  are  evaluated  from  their  physical  properties,  and  hardware  parameters 
of  the  device.48 
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(17) 


*«  («A>)  [(MsTsPp)  / (MpT„Ps)\m 

ftfun=(dAs)/  SDAP  (18) 

n - flow  rate,  M - molecular  weight  of  gas  mixture,  T - temperature,  P - pressure,  d - height  of  primary  flow  channel,  A 
- flow  cross-section,  D - orifice  diameter  of  secondary  flow  (index  p - primary  flow,  index  s - secondary  flow).  The 
best  mixing  conditions  resulting  in  the  maximum  output  power  were  at  the  values  of  between  0.8- 1 .0. 

State  of  the  art  of  the  performance  characteristics  of  the  COIL  device  operated  in  our  laboratory  are  summarized  in 
Table  3.  The  laser  operation  has  not  been  assessed  as  enough  efficient  yet.  The  reason  for  a rather  low  chemical 
efficiency  (measured  output  power)  is  most  probably  a bad  quality  of  the  resonator  mirrors  used  so  far  in  this  laser 
system. 

Tab.  3. 


Laser  power 

430  W 

Running  time 

1.5  min 

Chlorine  flowrate 

37.8  mmol/s 

Primary  helium  diluent 

80  mmol/s 

Generator  pressure 

8 kPa  (60  Torr) 

Chlorine  utilization 

0.97 

Starting  BHP  molarity 

6.7  M H02' 

02(1Ae)  measured  yield 

0.72 

Pressure  upstream  the  nozzle 

3.8  kPa  (28.5  Torr) 

Laser  cavity  pressure 

380  Pa  (2.8  Torr) 

1 2 /02  flow  rate  ratio 

0.029 

Mirrors  reflectivity 

0.9995/0.981 

Mode  length 

3.7  cm 

Small  signal  gain  (Rigrod) 

O.OBcnV1 

Saturation  intensity 

5.4  kW/cm2 

Chemical  efficiency 

0.12 

5.2.  Small  signal  gain  distribution  and  cavity  temperature  evaluation 

A spatial  distribution  of  the  small  signal  gain,  and  evaluation  of  cavity  temperature  was  investigated  by  means  of  the 
Iodine  Scan  Diode  Probe  diagnostics  developed  for  a COIL  at  the  PSI  , and  used  first  at  the  AFRL.50  The  line-shape  of 
the  COIL  gain  profile  was  scanned  in  frequency  by  a diode  laser  of  narrow  line  width  operating  in  the  region  of  the 
COIL  transition  frequency.  The  small  signal  gain,  g(v0),  was  determined  using  the  equation47 

I(v0)  = Io(v0)  exp(g(v0)  L)  (19) 

I^Vq)  - line-center  probe  laser  intensity  entering  the  gain  region,  I(v0)  - line-center  intensity  after  amplification  through 
the  gain  region,  L - gain  length. 

The  temperature  in  terms  of  the  Doppler  full  width  at  half-maximum,  Avd,  was  evaluated  from  the  relation47 

Avd  (7)  = 2 v0  / c {In  2 (2  kTIM)} 1/2  (20) 

v0  - central  frequency,  c - light  velocity  within  medium,  k - the  Boltzmann’s  constant,  M - mass  of  medium  species. 

Performed  mapping  of  the  gain  distribution  through  the  gain  region  was  aimed  at  the  inspection  of  mixing  conditions  in 
the  cavity  at  different  experimental  conditions,  and  bearing  out  our  guess  at  the  reason  of  rather  low  chemical  efficiency 
of  this  device  (mirrors  quality).  The  results  presented  below  were  gathered  for  the  primary  gas  flow  of  40  mmol  Cl2/s 
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and  80  mmol  Heprim/s,  the  secondary  gas  flow  of  40-80  Hesec  mmol/s  and  1-2  mmol  I2/s.  The  plenum  pressure  was  3.4  - 
4 kPa  (25-30  Torr),  the  cavity  pressure  measured  at  the  upper  wall  in  the  location  of  the  resonator  axis  was  300-390  Pa 
(2.2-2. 8 torr).  The  2.5  mm  i.d.-beam  of  the  probe  laser  was  passed  through  the  gain  medium  parallel  to  the  optical  axis. 
A two-pass  configuration  of  the  probe  beam  was  chosen  for  the  data  extraction.  The  probe  beam  emitter/detector  unit 
and  the  mirror  system  for  beam  alignment  were  mounted  on  an  assembly  of  the  motorized  linear  positioning  equipment 
that  allowed  fast  and  precise  movements  either  in  horizontal  direction  (x  axis)  along  the  gas  flow  or  vertical  direction  (y 
axis)  perpendicular  to  the  gas  flow.  The  grid  of  measuring  positions  is  presented  in  Fig.  11. 


Fig.  11.  Scheme  of  probe  beam  scanning  in  vertical  and  horizontal  position 

Duration  of  each  scan  was  13.3  seconds,  and  the  gain  data  were  recorded  within  60-100  seconds’  hot  runs.  One  example 
of  the  gain  distribution  curves  cross  the  gain  region  in  the  vertical  direction  measured  at  three  values  of  1 2 flow  rate  is 
illustrated  in  Fig.  12.  A parabolic  shape  of  the  gain  distribution  curve  indicates  that  the  gas  flow  velocity  across  the  gain 
region  was  inhomogeneous  resulting  in  a slower  gas  flow  at  boundary  layers  (0.5  to  1.5  mm  thick)  along  the  cavity 
walls.  In  these  layers,  a remarkable  quenching  of  excited  I atoms  occurred.  Molecular  iodine  may  be  also  present  in 
these  layers  which  may  consume  02(iAg)  in  its  dissociation.  The  maximum  gain  was  attained  at  the  iodine  flow  rate  of 
1.5  mmol  I2/s.  Similar  parabolic  distribution  of  gain  was  obtained  for  different  distance  from  the  nozzle,  but  the 
maximum  gain  values  were  attained  at  different  iodine  flow  rates  under  constant  other  conditions.  An  insufficient 
mixing  of  primary  and  secondary  flows  (inefficient  penetration)  could  be  easily  detected  by  this  diagnostics  (see  Fig. 
13).  The  gain  distribution  record  in  this  case  (40  mmol  Hesec/s  only)  showed  two  peaks  and  minimum  located  in  the  cavity 
center,  which  can  be  explained  by  a shortage  of  iodine  is  in  the  center  of  the  cavity. 


distance  from  center,  cm 


distance  from  center,  cm 


Fig.  12.  Gain  distribution  in  vertical  direction,  5.5  cm 
downstream  nozzle  throat  (resonator  axis);  for  three 
I2  flow  rates;  80  mmol  Hesec/s. 


Fig.  13.  Gain  distribution  in  vertical  direction,  5.5  cm 
downstream  nozzle  throat  (resonator  axis); 

1 .5  mmol  I2/s,  40  mmol  Hesec/s. 


A plot  of  measured  gain  on  evaluated  relative  penetration  factor  is  shown  in  Fig.  14.  The  small  signal  gain  of  ~l%/cm 
was  typical  for  the  supersonic  COIL  device  operated  in  our  laboratory  at  optimal  mixing  conditions.  Comparing  it,  for 
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example,  with  the  gain  attainable  on  the  RADICL  device  (-1,2  %/cm),  and  considering  also  twice  higher  cavity 
pressure  in  this  laser  system  (5.8  Torr),  there  is  a prospect  to  enhance  the  gain  of  our  COIL  driven  by  jet  SOG  when  the 
cavity  pressure  is  increased. 


Fig.  14.  Gain  in  the  cavity  center  in  relation  to  calculated 
relative  penetration  factor;  5.5  cm  from  the  nozzle 
throat  (resonator  axis),  1 .5  mmol  I2/s,  80  mmol  Hscc/s. 

The  gas  temperature  calculated  from  these  measurements  was  not  homogeneous  across  the  whole  cavity.  In  the 
boundary  layers  when  the  probe  beam  recorded  weak  emission  or  even  absorption,  the  gas  temperature  was  much 
higher  (>350  K)  than  in  the  cavity  center  (175  - 205  K).  A low  gas  velocity  and  high  rate  of  exothermic  processes  per 
unit  length  in  the  boundary  layers  near  the  walls  to  the  cavity  center  can  explain  it. 


6.  SUMMARY 

The  above-presented  results  have  been  extracted  from  a number  of  publications  and  reports  that  have  come  out  from  the 
COIL  Laboratory  in  Prague.  We  hope  that  our  current  and  future  research  on  COIL  related  problems  contributes  and 
will  contribute  to  the  advanced  COIL  technology. 
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